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Abstract 



Making use of an identity transformation independent of a nuclear model, we rep- 
resent the 2i//3/3-amplitude as a sum of two terms. One term accounts for most of 
the sensitivity of the original 2i//3/3-amplitude to g'^^ for realistic g'^^ ~ 1 (with g'^^ 
^ ' being the ratio of the triplet and singlet p-p interaction strengths) and is deter- 

Ch ■ mined by a specific energy-weighted sum rule. The sum rule depends only on the 

particle-particle residual interaction (being linear function of g'^^ in the QRPA) and 
passes through zero at the point g'^^ = 1 where the Wigner SU(4) symmetry is re- 
stored in the p-p sector of the Hamiltonian. The second term in the decomposition 
5^ , of the 2z^/3/3-amplitude is demonstrated within the QRPA to be a much smoother 

function for the realistic values of g'pp than the original 2i^/?/?-amplitude. This term 
is mainly determined by the intensity of the spin-orbit interaction of the nuclear 
mean field. Thus, the analysis of the present work reveals the reasons for the sen- 
sitivity of the 2z^/3/3-amplitude to different components of the nuclear Hamiltonian 
and thereby can help in constraining nuclear model uncertainties in calculations of 
the amplitude. 
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1 Introduction 



Numerous calculations of the 2i//3/3-amplitude within the QRPA have revealed 
its sensitive dependence on the strength Qpp of the charge-exchange triplet 
particle-particle (p-p) interaction (see, e.g., reviews [1]). The amplitude passes 
through zero at a realistic value of the interaction strength. However, the 
mechanism of such a dependence found for the first time in the pioneering 
work [2] almost 20 years ago still remains unexplained. The paper [2] already 
contained a clue that the vanishing of the amplitude may be related to the 
restoration of the Wigner spin-isospin SU(4) symmetry in nuclei. In an exactly 
solvable schematic model the authors of [2] also showed that this kind of 
dependence is not an artifact of the QRPA but shows up in the exact solution. 

There have been several semiquantitative attempts to use the concept of 
the dynamically broken SU(4) symmetry as a basis for describing the 21/(3(3- 
amplitude [3]. Since the amplitude is zero when the SU(4) symmetry is exact, 
it seems natural to express the amplitude explicitly in terms of those parts 
of the nuclear Hamiltonian H which are responsible for the violation of the 
symmetry. To achieve this goal, a model-independent transformation of the 
amplitude was suggested in [4]. Unfortunately, the quantitative realization of 
this idea in [4] was based on a simplified model of independent quasiparti- 
cles and, therefore, could not address the question of the sensitivity of the 
2i//9/?-amplitude calculated within the QRPA. That was the reason why the 
work [4] faced strong criticism claiming that "the formalism fails in describing 
the double beta decay observables" [5] . 

To refute the criticism of [5], in present work we substantially develop the 
approach [4] by incorporating the QRPA into it. As in [4], the starting point 
for the analysis is model-independent, identity transformation of the 2v(3(3- 
amplitude: M21, — + S/cUg. That allows us to separate out a part pro- 
portional to an specific energy-weighted sum rule S depending only on the 
particle-particle residual interaction. This part is shown within the QRPA to 
account for most of the sensitivity of the original 2z//9/9-amplitude to g'pp for 
realistic g'pp ~ 1 (with g'^^ being the ratio of the triplet and singlet p-p inter- 
action strengths). The sum rule S can be expressed in the form that avoids 
the intermediate state summation and involves only the wave functions of the 
initial and final ground states. An analytical representation for S that can be 
given within the QRPA using a simple separable p-p interaction reveals its 
simple linear dependence on g'^^ with passing through zero at g'^^ = 1 where 
the SU(4) symmetry is restored in the particle-particle sector of the Hamilto- 
nian. The second part in the decomposition of the 2z//3/3-amplitude can 
be explicitly expressed in terms of those components of the nuclear Hamilto- 
nian which violate the SU(4) symmetry. As shown by our QRPA calculations, 
this part behaves much smoother for the realistic values of g' and is mainly 
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determined by the spin-orbit part of the nuclear mean field as was conjectured 
by Rumyantsev and Urin [4]. 

Thus, we think that our present results contribute to understanding of general 
physical properties of the 2i//3/3-decay amplitude and bring new opportunity 
to constrain nuclear model uncertainties in calculations of the amplitude. In 
particular, it is shown that the sensitivity of the amplitude to the p-p inter- 
action strength has its natural physical origin in the violation of the SU(4) 
symmetry in the particle-particle sector of a realistic nuclear Hamiltonian and 
has nothing to do with the shortcomings of the QRPA itself. Also an impor- 
tant influence of the spin-orbit mean field on the calculation results that has 
received rather little attention in the preceding publications is emphasized. 
Besides, the sum rule S can be easily calculated in different nuclear models, 
for instance, in shell model, thereby providing a test of the quahty of the 
QRPA calcuations. 

Note, that we do not claim (as well as the authors of [4] did not) any model- 
independent calculations of the 2 Z//9/9- amplitude (that was misunderstood by 
authors of [5]) rather we are dealing with model-independent transformation of 
the 2i//3/3-amplitude. Nevertheless, as shown in the paper this identity transfor- 
mation allows us to shed some light on the general properties of the amplitude 
calculated within the QRPA. 

The paper is organized as follows. In section 2 an identity transformation of 
the 2i//3/3-amplitude which is valid independently of a nuclear structure model 
is performed. The transformation turns out to be a power tool for the analysis 
of the 2z//3/5-dccay amplitude because it allows one to incorporate explicitly 
the symmetry properties of the nuclear Hamiltonian. The model Hamiltonian 
consisting of a mean field, the Landau-Migdal particle-hole interaction and 
a separable particle-particle interaction is specified in section 3. The identity 
transformation separates out a part proportional to a specific energy-weighted 
sum rule for which an analytical representation is given within the QRPA 
in section 4. Calculation results for ^^Ge, ^°°Mo and ^^°Te are presented in 
section 5 and we conclude in section 6. 



2 Model-independent transformation of the 2i//3/3-ampIitude 

We adopt here the same fine of reasoning as presented in [4] . To make the fol- 
lowing derivation of the new representation for the 2z//9/9-decay amplitude 
more transparent, we will use the same approach for the Fermi and the 
Gamow- Teller transitions. Let G^"^^ = Y,a9^^ be the operators of allowed 
Fermi (y(±) = r^^)) and GT (^(±) = crr^^^) (3 - transitions. 
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We start our derivation with writing down the usual expression for the 2u(3 (3 
decay amplitude (see, e.g., reviews [1]): 

M2. = ^^. (1 



Here, Qg — {f\&~^\s), fs = {s\&^^\i) are the one-step matrix elements of 
the operator G^"^ between the intermediate states |s) of an isobaric nucleus 
{N — 1, Z + 1) and the ground states \i) and |/) of a parent {Z,N) and 
the final {Z + 2,N — 2) nucleus, respectively, with Es, Ei and Ef being the 
corresponding energies and uJs = Es — {Ei + Ef)/2 being the excitation energy 
of an intermediate state measured from the mean value of the ground-state 
energies. 

Then wc perform the following identical splitting of the amplitude (1) into 
two parts: 



^9 



s 



Y^^sOsfs-, 



(2) 

(3) 

(4) 



where Ug is an arbitrary energy that was identified with the energy uiq of the 
corresponding giant resonance in [4]. It should be so in the hmit of exact SU(4) 
symmetry, but for the realistic situation one can only say that uig should be 
rather close to loq- Note, that the second term in (2) is proportional to S 
having the form (4) of a specific energy-weighted sum rule. Bearing in mind 
that 



Js = {s\ [h,G'--^] \i)-{Ef-E,)fs/2, 
Qs^ifl \h,G^-A \s) - {Ef - E,)gs/2 , 



(5) 



where H is the nuclear Hamiltonian, one can get the following expression for 
S: 



S^l{f\[G^-\[H,G^~y 



(6) 



Here we do not support the claim of [4] that S should vanish (see Eq. (5) of the 
citation). That would be true only if the p-p interactions in the Hamiltonian 
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could be neglected. But it has been known for a long time that the p-p inter- 
actions play an important role in the nuclei far from closed shells, especially 
for describing 2i//3/3-decay amplitudes [2]. Actually, the authors of [5] provided 
direct evidence by their QRPA calculations that S does not vanish (except for 
one point in gpp dependence, see Fig. 1 of the reference) . 

We will show in Section 4 that one can get within the QRPA an analytical 
expression for S depending linearly on Qpp. The point where S = then corre- 
sponds to the restoration of the SU(4)-symmetry in the particle-particle sector 
of the Hamiltonian (i.e. commutes with the p-p part of the Hamiltonian 
for this gpp). 

To conclude this section, we show, following [4], that M!^,^ can be written as 



gsfs 



00. 



2 ^ 



(7) 



where/, = (s|T4"^|i) = -{oug - u;s){s\G^-^t) and^, = {f\Vr\s) = -{oug + u;s){f\G^-^\s 
are the matrix elements of the operator V^~^ defined as: 



Vq — 



Ef — Ei 



(8) 



It is noteworthy that M2j^ acquires contributions from all terms of the Hamil- 
tonian violating the SU(4) symmetry (in contrast to the sum rule S (6) which 
is determined only by the residual p-p interaction). On the mean field level the 
most important source of the violation is the spin-orbit mean field as was con- 
jectured by Rumyantsev and Urin [4]. We will see in Section 5 how strongly 
is affected by varying the intensity of the spin-orbit mean field in the 
calculations. 

Thus, the final expressions (6) and (7), (8) contain explicit information about 
symmetry properties of the nuclear Hamiltonian in terms of the corresponding 
commutators. Also we would like to stress again that the results of this section 
— the transformation of the 2z//3/?-amplitudc Eqs.(2)-(4) and the representa- 
tions for S and Eqs.(6)-(7) — are quite general and do not rely on any 
nuclear model. 



3 Model Hamiltonian 



For a quantitative realization of the relationships of the preceding section one 
has to specify the nuclear Hamiltonian. We adopt here the realistic model 
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Hamiltonian of [6] consisting of a mean field U{x). the Landau-Migdal zero- 
range forces for describing the particle-hole interaction Hp-^ and a separable 
particle-particle interaction in both the neutral and the charge-exchange chan- 
nels Hp-p. Only isoscalar (central Uo{r) plus spin-orbit Uso{T)crl) part of the 
mean field U (x) is treated phenomenologically while the isovector (symmetry 
potential) and the Coulomb parts are calculated consistently in the Hartree 
approximation (see [6]). The explicit expression for the Landau-Migdal forces 
in the charge-exchange channel is [7]: 

Hp.n = 2 ^(Fo + F,(Ta<T,)rt^Tt^S{r, - n) + h.c. (9) 

a>b 

where the intensities of the non-spin-flip and spin-flip parts of this interaction, 
Fq and Fi respectively, are the phenomenological parameters. We choose the 
p-p interaction in the following separable form (as used in the BCS model) in 
both the neutral (pairing) and charge-exchange particle-particle channels: 

J/J, j3=p,n 
P'=p,n 

where Gj are the intensities of the particle-particle interaction. Here we con- 
sider only the nucleon pairs with total angular momentum and parity J'^ = 
O"*", l"*" {L = 0,S — J). The annihilation operator for a nucleon pair P'^'* is 
given as follows: 

P^^' = T.UU'x')P^^',xx' ; P^^'M' - 4a)(a') T^iMj^xruj^yrny^yrn'^rn , (11) 

AA' mm' 

where axmiaxm) annihilation (creation) operator of a nucleon in the 

state with the quantum numbers Am (m is the projection of the particle 
angular momentum) and the single-particle radial wave functions r'^xxi^)^ 
ix^rXu) ^ IXn{r)xu{r)dr, = \/2i7r + l^WM, = 75((7r)||(7||(z/)). 

Note that, apart from the p-p interaction, the model Hamiltonian has proven 
its reliability by successful continuum- RPA calculations of many structure and 
decay properties of charge-exchange giant resonances (see, e.g., [8]). 

We use the Bogolyubov transformation to describe the nucleon pairing in the 
neutral channels in terms of quasiparticle creation (annihilation) operators 
'^Am('^Am) (see, e.g., [9]). As a result, the following model Hamiltonian is ob- 
tained to describe Fermi and Gamow-Teller excitations within the pn-QRPA: 

H = Ho + Hp-H + Hp-p + i^r.N + l2pZ; = E E^^(«AJ+«Am (12) 

P=p,n Am 
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Here, — yi^xY + quasiparticle energy, — — Hfj^ and 

are the chemical potential and the energy gap, respectively, which are 
determined from the BCS-model equations: 

^/3 = E(2jf + l)Kf; A;3 = GoE(2jf + l)«K (13) 

A A 

with vl = i(l - ^) and ul = l- vl 

The explicit expressions for the p-h and p-p interactions in terms of the quasi- 
particle (pn)-pair creation and annihilation operators are: 



Hp-h = E ^ E (x.x.x.'X.')(Qi'.rQX', (14) 

J/i nun'u' 

Hp-p = -E^^ E ix.Xu){x.'X.'){PiTPX', (15) 

J/Lt TTI/Tt'i/' 

where A-^i^ = E {Jf^\jnmjum')a^rn'a7rm, Kv = {-^y^'^Kv'^ ix^XuXn'Xw) = 

mm' 

I XttXi^Xtt'Xi/'^^^c?^- Note, that the radial integrals in the representation (15) 
for the separable p-p interaction should be kept rather than put equal to 1 
as in many applications of the separable forces. This is due to a slight non- 
orthogonality of proton and neutron wave functions having the same quantum 
numbers. 

Our model Hamiltonian preserves the isospin symmetry provided that a self- 
consistency condition relating the isovector part of the mean field (symmetry 
potential) to the Landau-Migdal parameter Fq is fulfilled (see [6], the p-p in- 
teraction (10) is isospin invariant by construction). That was proved by the 
direct continuum-QRPA calculations of the properties of the isobaric analog 
states in [6] . Note, that the usual discretized version of the pn-QRPA which 
we are going to exploit here violates a bit the isospin symmetry due to a re- 
stricted basis of the particle-hole configurations. Apart from this source of the 
violation, we shall introduce another, usual one, namely, we allow the pair- 
ing constants for proton and neutron subsystems (Gop and Gon, respectively) 
to have different values in order to reproduce the experimental pairing gaps. 
These approximations certainly affect the quality of the calculated Fermi ma- 
trix elements that is worth to discuss elsewhere. However, they seem to be 
of rather minor importance for the description of the Gamow- Teller matrix 
elements that we are pursuing here. 
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4 Analytical representation for S 



In this section we show that a simple and elegant analytical representation for 
the EWSR S (4) introduced in Section 2 can be obtained using the Hamilto- 
nian of the preceding section. To do that, one has to calculate the expectation 
value of the double commutator (6). The equation of motion for the GT op- 
erator can be derived within the QRPA [4,6]: 



{01-% = 



Gon — Gi 



On 



'PP 

'^P-^ pn ■ 



'Op_ 



(16) 

(17) 



Op 



The rhs of (16) is the sum of particle-hole M and particle-particle 



ph 



Ul M operators originating from the commutator of & ^ with the p-h 



pp 



part (including the mean field) and the p-p part of the Hamiltonian, respec- 
tively. The explicit representation for (f^^ ^ can be found in [4,6] but it is 

= Odue 



ph 



irrelevant for the current derivation of S since (/[ 
to T^~^^ = 0. Hence, only p-p part of the Hamiltonian contribute to -S": 

A„Ap Gqji + Gqp 



(/I 



V / pp 



3(1 - gl 



pp' 



GonGop 



(18) 



where we have introduced the parameter g'^^ — 2Gi/{Gon + Gop) describing 
the relative strength of p-p interaction in the triplet spin channel with respect 
to the singlet (pairing) one (not to be confused with the parameter Qpp used 
in the calculations with realistic nucleon-nucleon interactions to renormalize 
the p-p part of them). In getting (18) we have used an approximation that the 
BCS vacua are the same for the initial and final ground states. For the Fermi 
transitions factor 3 in (18) should be omitted and g'pp — 2Go/{Gon + Gop). 

In the case of the restored isospin symmetry in the p-p channel, i.e. Gon — 
Gqp = Go, the expression (18) can be simplified further: 

S-S{l-g'J^. (19) 



The linear dependence of the sum rule S on g'^^ with crossing zero at the 
point of the restoration of SU(4)-symmetry {g'^^ = 1) in the p-p channel is 
of great importance. It shows analytically that at least part of the 2z//3/3- 
decay amplitude reveals a quite general sensitivity to g' which is not really 
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restricted to the used approximations like the QRPA and BCS. In fact, any 
reahstic nuclear Hamiltonian should contain singlet and triplet p-p interactions 
and any model of the ground states of open-shell nuclei should account for 
the pairing effects, in particular, the increased probability to pick up a pair of 
nucleons in the O"*" state which is described by A in the BCS model. Therefore, 
we believe that S will be linearly decreasing with respect to g'^p in any realistic 
model (including modern shell models) while the slope of the dependence and 
its zero may be slightly different from those given by (18), (19). ^ 



5 Calculation results 

To illustrate the method suggested, we present in this section some calcu- 
lation results obtained within the pn-QRPA for ^^Ge, ^''^Mo and ^^''Te. Pa- 
rameterization of the isoscalar part of the mean field Uo{x) along with the 
values of the parameters used in the calculations have been described in de- 
tails in [10]. Here we increase the depth of the central part Uo{r) by about 
10 MeV to get a large enough basis of the bound single-particle states (this 
change of a isoscalar parameter should not afi:ect much the isovector quan- 
tities like the Gamow- Teller matrix elements). The dimensionless intensities 
fj of the Landau- Migdal forces of Eq. (9) are chosen as usual: Fj — fj ■ 300 
MeVfm'^. The value fo = f' = 1.0 of the parameter /' determining the symme- 
try potential according to the selfconsistency condition is also taken from [10], 
where the experimental nucleon separation energies have been satisfactorily 
described for closed-shell subsystems for a number of nuclei. The dimension- 
less intensity /i = g' of the spin-isovector part of the Landau- Migdal forces is 
chosen as g' = 0.6 that allowes one to reproduce rather well the experimental 
energies of GTR in the nuclei in question. The value g' = 1.0 would correspond 
to the restoration of the SU(4) symmetry for p-h interaction in our model. 

The strengths of the monopole pairing particle-particle interaction Gon, 
are chosen to reproduce the experimental pairing energies in the considered 
nuclei. The BCS basis has 9 levels (oscillator shells A^=3,4) for ^^Ge, 10 levels 
(oscillator shells A^=3,4 plus the l/iii/2 orbit from N = 5) for ^°°Mo, and 12 
levels (oscillator shells A^=3,4 plus the 1/^9/2,11/2 + 2/7/2 orbits from N — 5) for 
130 Pqj. g^j.jj basis systems the reahstic values are Gon ~ Gop 12 MeV /A 
with at most a variation of 20%. Hereafter, we will refer to this set of the 
pairing constants as set I. We have checked that increasing the BCS basis by 
adding iV=l,2 shells change the results only slightly. Therefore, in addition, 
we choose artificially large Gon = Gop — IS MeV /A (referred to as set H) as a 
test case in order to emphasize some important points of our approach. 

^ This our conjecture has been confirmed by direct shell-model calculations of an 
anonymous Referee to whom we are indebted. 
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The strength of the triplet spin-spin particle-particle interaction Gi is con- 
sidered as a free parameter and aU dependencies are plotted as functions of 
the parameter g'^^ = 2Gi/(Gon + G^p). Note, that the basis for the QRPA 
diagonalization consists of all the bound s.p. states, thereby containing all the 
main transitions among the spin-orbit doublets with the same radial quantum 
numbers. The explicit representation for the set of the QRPA equations can 
be found elsewhere (see, e.g., [1]). Note, that both j3~ and branches to 
construct the 2i//?/3-amplitude are calculated within the QRPA in the same, 
initial for the 2i//3/9-decay, nucleus, so we adopt here the same approximation 
as in [2]. 

In Fig. 1 we show S{g'p^) calculated according to the definition Eq. (4) (solid 
line) and its analytical representation Eq. (18) (dashed line). The horizontal 
dotted line corresponds to 5" = 0. Also the dependences corresponding to the 
pairing parameter sets I and II are plotted, where the latter is steeper due 
to the larger Ap, A„. In all cases the deviation of the numerical results from 
the analytical ones (with the largest deviation observed in ^^'^Mo) can barely 
be seen despite the fact that our QRPA bases do not contain the whole p- 
h spectrum as assumed in deriving Eq. (18). It gives one more confirmation 
that the QRPA is capable of reproducing well various sum rules. Larger pairing 
correlations produce larger deviations due to a larger contribution of the p- 
h configurations discarded in the truncated basis. It is noteworthy that the 
calculations of S in ™Ge using the G-matrix generated residual interaction 
have also revealed the linear dependence S{gpp) with almost the same value 
S{gpp — 0) (see Fig. 1 of [5]). It can also be seen from Fig. 1 that the realistic 
dependence S{g'p^ is almost the same for different nuclei. It can be understood 
from the phenomenological proportionalities A oc A~^/'^ and Go oc put 
into S cx A„Ap/Go- 

In Fig. 2 we plot M2u{g'pp) (solid hue) calculated according to Eq. (1), S/u'^g 
(dashed line) and their difference M'^i^ig'pp) (dot-dashed fine). The value of Ug 
is calculated as the mean energy of the Gamow- Teller strength distribution in 
the /9^-channel. For each nucleus we have two sets of calculations - left and 
right panels contain the results corresponding to the pairing constant sets 1 
and II, respectively. One can see from the figure that the dependence M2^{g'pp) 
is smoother than the original one, M2v{g'pp) (at least for not-too-large values of 
g'pp where the dependencies become essentially non-linear in the parameter). 
Furthermore, by going to the unphysical set 11 of the pairing parameters, one 
can see that the main change in the amplitude of M2u{g'pp) can be attributed 
to the change of S{g'pp) with a smaller change of M2i,{g'pp). As for the function 
S/ujg, it shows basically a hnear dependence on g'^^ controlled by S. 

Finally, we have tested the main conjecture of [4] that M21, (7) is mainly 
determined by the spin-orbit mean field contribution to the commutator (8). 
We plot the corresponding results in Fig. 3. The sohd and dash-dotted lines 
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represent Ml^^lg'pp) calculated according to the definition (3) for 100% and 50% 
of the realistic strength of the spin-orbit interaction, respectively. The dashed 
and dotted lines represent Mi^j^^g'^p) calculated according to the (7), (8), taking 
only spin-orbit contribution to the latter equation (see [4] for the details), 
also for 100% and 50% of the realistic strength of the spin-orbit interaction, 
respectively. One can see from Fig. 3 that M2^{gpp) is indeed very sensitive 
to the change of Uso- The approximation of M2^{g'pp) by (7), (8) is fairly good 
for the realistic g'^p in the vicinity of unity (a recent discussion of the realistic 
gpp can be found in [11]). Thus, in that region Ml^^^g'^p) should scale almost 
quadratically with the magnitude of Ugo and the realistic choice of the spin- 
orbit interaction intensity is very important. It is also noteworthy that the 
schematic model QRPA calculations [2] revealed almost linear dependence of 
the 2i//3/3-amplitude on g'^^ that is quite understandable because the authors 
neglected the spin-orbit splitting in the model ^ . 



6 Conclusions 

In this paper, using an identity transformation, we have represented the 2u(3(3- 
amplitude as a sum of two terms. One term has been shown to be proportional 
to an specific energy- weighted sum rule depending only on the particle-particle 
part of the Hamiltonian. This part has been shown within the QRPA to ac- 
count for most of the sensitivity of the original 2z//5/5-amplitude to g'^^ for 
realistic g'^^ ~ 1. The sum rule has been found within the QRPA to depend 
linearly on the ratio g'pp and to vanish at the point g'pp = 1 where the Wigner 
SU(4) symmetry is restored in the p-p sector of the Hamiltonian. The ana- 
lytical results for the sum rule have been reproduced very well by the QRPA 
calculations using truncated p-h bases. Thus, we have shown that the sen- 
sitivity of the amplitude to g'pp is due to the pronounced effect of violation 
of the SU(4) symmetry in the p-p sector of a reahstic nuclear Hamiltonian. 
The second term, M!^^, in the representation of the 2i//3/3-amphtude has been 
shown within the QRPA to be a much smoother function for realistic values 
of g'pp than the original 2z//5/5-amplitude. We have confirmed the conjecture of 
Rumyantsev and Urin [4] that is mainly determined by the intensity of 
the spin-orbit interaction of the nuclear mean field and scales approximately 
as square of the intensity for the reahstic values of g'pp. Therefore, the real- 
istic choice of the spin-orbit mean field is very important for describing the 
2z//?/5-amplitude. Finally, the analysis of the present work has revealed the 
reasons for the sensitivity of the 2i//3/3-amplitude to different components of 
the nuclear Hamiltonian and thereby can help in constraining nuclear model 
uncertainties in calculations of the amplitude. 

^ unfortunately, it is not possible to compare our results with the values of the 
2i^/3/3-amplitude from Fig. 2 of [2] because the ordinate was given in arbitrary units. 
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Fig. 1. The energy- weighted sum rule S calculated according to the definition, 
Eq. (4), (solid line) and its analytical representation Eq. (18), (dashed line), for 
the pairing parameter sets I and II (the latter is steeper). 
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Fig. 2. The decomposition of the calculated M2u{g'pp) (solid line) as the sum of 
S/LOg (dashed line) and M2y{g'pp) (dot-dashed line). The two sets of calculations for 
each nucleus are contained in the left and right panels corresponding to the pairing 
constant sets I and II, respectively. 
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Fig. 3. M2j^ calculated according to the definition (3) for 100% (solid line) and 50% 
(dash-dotted line) of the realistic strength of the spin-orbit interaction. The dashed 
and dotted lines represent calculated according to the (7), (8), taking only the 
spin-orbit contribution to the latter equation for 100% and 50% of the strength of 
the spin-orbit interaction, respectively. 
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